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ABSTRACT 

The magnetic fields of white dwarfs distort their shape generating an anisotropic moment of 
inertia. A magnetized white dwarf which rotates obliquely relative to the symmetry axis has a 
mass quadrupole moment which varies in time, so it will emit gravitational radiation. LISA may 
be able to detect the gravitational waves from two nearby, quickly rotating white dwarfs. 



Introduction 



Compact binaries consisting of a white dwarf and a late main-sequence star such as AM Her and DQ Her 
cataclysmic variables (CVs) are one of the primary sources of gravitational radiation observable by the LISA 



mission (Bender et al. 1998). As the two stars orbit each other with a period of a few hours, the energy 
of their orbit is slowly lost to outgoing gravitational waves. Additionally, the white dwarfs themselves are 
rotating. For AM Her binaries, the white dwarf rotation period is approximately synchronized with the 
orbital period. The white dwarf in a DQ Her binary rotates ten or more times faster than the binary orbit. 
In particular, the DQ-Her, white dwarfs AE Aquarii and WZ Sagittae rotate with a periods of 33 and 28 
seconds respectively. If these two stars have a significant quadrupole moment that is not aligned with the 
rotation axis, they will be sources of observable gravitational radiation. 

Several authors have examined whether quickly rotating white dwarfs on their own will be important 



sources of gravitational waves. Lindblom (1999) explored the possibility that gravitational waves can drive 

; [Friedman (19831 



an instability in the r-modes of near-critical rotating white dwarfs analogous to process that 
Cutler fc Lindblom (19"87)t [Owen et al. (1998)| , [Lindblom, Owen fc Morsink (1998)t [Andersson, Kokkota; 
fe Schutz (1999) have argued may operate in rapidly rotating neutron stars. Lindblom found that this 



instability would take several gigayears to grow in the observed DQ Her objects. Furthermore, the r-mode 
instability can only be excited in hot (T > 10® K), massive white dwarfs (M > 0.9 Mq) rotating within 
twenty percent of their critical period; therefore, it is unlikely that this instability plays an important role 
in white dwarfs observable by LISA. 

The magnetic field of neutron stars may induce a sufficiently large quadrupole moment to affect their 
spin-down ( Shapiro fc Teukolsky 1983| ) and to produce observable gravitational radiation (Bonazzola & 
Gourgoulhon 199(:, [Brady et al. 1998 ). The internal magnetic fields of white dwarfs distort their structure 
similarly {e.g. Ostriker fc Hartwick 1968| ). If a quickly rotating white dwarf harbors a strong internal magnetic 
field, it may be an important source of gravitational radiation. 



The fastest rotating white dwarfs have periods of several tens of seconds while typical white dwarfs rotate 
in about hour. This places them in the most sensitive region of the spectrum for the LISA interferometer. 
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2. Magnetically Induced Distortion 



The magnetic field supplies an anisotropic pressure. The tensor virial theorem quantifies the effect of 
the field on the structure of the star. If the field points along the z-axis and the total magnetic energy (M) 
is much less than the internal energy (H), 
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where Wij is the potential energy tensor (c./. Binney & Tremaine 1987), W — Wu is the total potential 



energy of the star, and z is the symmetry axis of the magnetic configuration 

The lowest order term in a multipole expansion of the gravitational radiation in the wave zone 
{i.e. d ^ c/il) of a slowly moving source is the quadrupole term. Because general relativity only modifies 
the structure of white dwarfs slightly, the reduced quadrupole moment (/y ) is well approximated by the 



Newtonian limit ( Misner, Thornc fc Wheeler 1973 ) 
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If one assumes that the internal field is uniform and that the isodensity surfaces of the distorted star are 



similar ellipsoids ([Binney fc Tremaine 1987| ), this yields 
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where Iq is the moment of inertia of the unperturbed spherical star. A negative value of Izz indicates that 
the configuration is oblate. 

I I 

Chandrasekhar & Fermi (1953) treated a similar problem by minimizing the total energy of a 

uniform density sphere with a uniform internal field and a dipolar external field by perturbing its surface 
axisymmetrically. They obtained 
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for small quadrupolar deformations. 



Detailed modeling indicates that the distortion also depends on the structure of the internal magnetic 
field. Specifically stars with a larger fraction of toroidal fields tend to be prolate while stars whose fields are 
dominated by a poloidal component are oblate. The distorted structure of the star is calculated by assuming 
that the equation of state is strictly barotropic, that both the magnetic field and rotation only slightly 
perturb the structure of the star and that the current and the field vanish at the surface ( Woltjer 196C| , 



Wentzcl 1961, Ostriker & Hartwick 1968). The conditions specify the field through an eigenvalue equation. 



The Chandrasekhar (1939) equation of state (with fj,e — 2) determines the unperturbed structure of the 
star. This equation of state provides an adequate approximation for this purpose and is relatively simple 



analytically compared to equations of state which include Coulomb and other corrections {e.g. Hamada & 
Salpeter 196l| ). 



Figure |] shows how the reduced quadrupole moment depends on the structure of the field and the 
central density of the white dwarf. As the toroidal field increases relative to the poloidal field, the star 
becomes more prolate. Additionally, as the central density increases the white dwarf becomes more centrally 
concentrated, and the magnetic field is less effective at distorting the star. 
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Independent of the configuration of the internal magnetic field, the reduced quadrupole moment is given 
approximately by 

4. « 0.15^/0 (5) 

for allowed eigenstates of the magnetic configuration. The moment of inertia of the unperturbed star (Jq) 
decreases from Q.2MR^ to O.OSMi?^ as the model becomes more centrally concentrated. 

The value of I^z differs in sign from the estimates (Equation ^ and ^) . The surface boundary condition 
forces the field to have a large toroidal component that tends to make the star prolate. The results which do 
satisfy the boundary condition may be extrapolated linearly to examine fully poloidal configurations (which 
have surface currents). A fully poloidal configuration has 



while a fully toroidal field has 



- 0.25^/0. (7) 



These distortions fall short of the simple estimates by an order of magnitude because the magnetic field is 
centrally concentrated, weakening its effect on the outer layers which dominate the moment of inertia. 

If the magnetic field is not aligned with the rotation axis of the star, the star will emit gravitational 
radiation. If the symmetry axis of the field makes an angle a with the rotation axis, the components of the 
gravitational wave at the Earth are 

/i+ = ft,o sina I - cos a sini cosi cos [f2(t — to)] — sinck cos^(2z) cos [2r2(i — io)] I (8) 

/ix = /lo sina — cos a sini cosi sin [r2(i — ig)] ^ sina cos i sin [2f2(< — io)]|' (9) 

= -^4.^ (10) 

where the inclination between the rotation axis and the line of sight is i, and d is the distance from the star 
to the detector (Bonazzola & Gourgoulhon 1996). 



Potential Sources 



White dwarfs with strong internal fields may have no surface field. However, it is more likely that strong 
internal magnetic fields lurk inside white dwarfs with large external fields. AM Her and DQ Her CVs consist 
of a magnetic white dwarf accreting from a late main-sequence companion. In AM Her stars, the accreted 
material flows directly onto the magnetic field lines of the degenerate star {e.g. Cropper 1990| ). The DQ Her 
stars are the white-dwarf analogue to low-mass X-ray binaries in which the accreted material forms a disk 



around the primary before flowing into the white dwarf magnetosphere {e.g. Patterson 1994). The white 



dwarfs in these two types of binaries rotate with periods ranging from tens of seconds to several hours. The 
strength of the surface magnetic fleld (for those stars where the magnetic field has been determined) ranges 
from 10-100 MG. 
I 1 



Patterson (1998) and Lasota, Kuulkers & Charles (1999) argue that the dwarf nova WZ Sge is indeed 



a DQ Her star during quiescence. As it is thought to be relatively nearby (Meyer-Hofmeister, Meyer & 
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Liu 1998), and it is the fastest rotating putative member of the DQ Her class of magnetic CVs, it may be 
an important source of gravitational radiation. 

About two percent of non-accreting white dwarfs are strongly magnetized with surface magnetic fields 



ranging from several tens of kilogauss to one gigagauss (Schmidt & Smith 1995, Putney 199(;, Anselowitz 



et al. 1999). The periods of these objects range from a dozen minutes to over a century, making them a 
disparate class. The more strongly magnetized and quickly rotating members of this group may be important 
gravitational wave sources. 

Table |^ lists the potential sources along with an estimate of the strength of the gravitational radiation 
at the Earth assuming that the white dwarfs have masses of O.6M0, the their RMS internal field is 10^^ G 
and a = 7r/2. This yields values oiM/\W\ = 4 x lO'^ and 4^, = 6 x IQ-^h = 1 x IQ-*^ g cni^. The isolated 
(non-accreting) white dwarfs are denoted by IWD. 

Figure ^ compares the expected gravitational wave amplitudes for these sources with recent estimates 
for LISA of instrument thresholds for achieving a signal-to-noise ratio of unity and five over a one-year 



observation (Bender et al. 1998|). Also, an estimate of the binary confusion noise and the strength of the 



radiation from the binary WZ Sge are depicted. Unless the internal fields of most of the potential sources 
are much larger than 10^^ G, or they are much less massive that O.6M0, they will not be detected by LISA. 
However, the two fastest rotating white dwarfs in the sample, WZ Sge and AE Aqr lie above the signal-to- 
noise unity line, so these objects may be detectable after a single year's integration. In fact, WZ Sge may 
appear at two frequencies corresponding to its orbital and rotational periods. 



Discussion 



The internal magnetic fields of white dwarfs are difficult to constrain observationally or theoretically. 
Since the magnetic field may vanish at the surface, it may leave no obvious trace in the observed spectra 
of the star and only affect the global properties of the star, such as its moment of inertia. The frequencies 
of the pulsation modes of the star may also depend on the structure and strength of the magnetic field 
(Jones et al. 198£, Lou 1995| ). However, asteroseismological kernels for white dwarfs tend to have little 
amplitude in the central regions of the star. If a mode has a node within a composition transition region, 
it is pinched between the interface and the surface, and its amplitude and frequency depends little on the 



properties of the material below the interface (Kawaler, Novikov & Srinivasan 1997). Furthermore, unless 



the white dwarf is hot, the interior is nearly degenerate, and the Brunt- Vaisala frequency almost vanishes. 



so the observed g- modes are damped ( Hansen &: Kawaler 1994 , Kawaler, Sekii fc Gough 199E ) through in 



the interior. These effects make it difficult to constrain the magnetic field strength in white dwarf interiors 
through asteroseismology. 

The strong fields of magnetic white dwarfs may be remnant of the magnetic fields of their progenitors 
on the main sequence or produced subsequently. As the core of the star forms a white dwarf, if the magnetic 



flux is conserved, the resulting remnant will have a similar ratio of magnetic to gravitational energy ( Ostriker 



fe Hartwick 1968); therefore, if the fields are relic, one can estimate their strength by studying stars along 



the main sequence. 

In main sequence stars it is difficult to exclude the presence of magnetic fields with a small fraction 
of the binding energy of the star. Helioseismological measurements of the sun place an upper limit on its 



internal magnetic field (Goode et al. 1996). The magnetic pressure throughout the solar interior is less than 
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Table 1: Potential White Dwarf Sources of Gravitational Radiation: Nearby magnetized white dwarfs with 



rotational periods less than 3.5'*. The data on the various types of sources was culled from: Cropper (1990) 



for AM Her objects, [Patterson (1994)|for DQ Her objects and [ Anselowitz et al. (1999)| for IWD. The data for 



RE J0317+853 are from Burleigh, Jordan & Schweizer (1999) and for WZ Sge are from Meycr-Hofmeister 



Meyer fc Liu (1998) The distances to AE Aqr and V471 Tau are from Ferryman et al. (1997) 



Source 


Period [s] 


Distance [pc] 


Bs [MG] 






h 




Type 


WZ Sge 


28 


50 


- 


4 


X 


10" 


- 


Dwarf Nova 


AE Aqr 


33 


102 


50 


1 


X 


10" 


01 

zo 


DQ Her 


V533 Her 


64 


1000 




4 


X 


10" 




DQ Her 


DQ Her 


142 


420 




2 


X 


10" 


zo 


DQ Her 


H 0253+193 


206 


200 




2 


X 


10" 




DQ Her 


GK Per 


351 


490 




3 


X 


10" 


26 


DQ Her 


YY Dra 


529 


155 




4 


X 


10" 


zo 


DQ Her 


V471 Tau 


555 


46 


— 


1 


X 


10" 


25 


DQ Her 


RE J0317+853 


726 


35 


660 


1 


X 


10" 


25 


IWD 


V1223 Sgr 


746 


400 




8 


X 


10" 


27 


DQ Her 


AO Psc 


805 


250 




1 


X 


10" 


26 


DQ Her 


RE 0751+144 


834 


400 




6 


X 


10" 


27 


DQ Her 


BG CMi 


913 


500 




4 


X 


10" 


27 


DQ Her 


FO Aqr 


1254 


300 




4 


X 


10" 


27 


DQ Her 


TV Col 


1910 


400 




1 


X 


10" 


27 


DQ Her 


TX Col 


1911 


500 




1 


X 


10" 


27 


DQ Her 


VZ Pyx 


2918 


250 




9 


X 


10" 


28 


DQ Her 


V1062 Tau 


3726 


1100 




1 


X 


10" 


28 


DQ Her 


EX Hya 


4022 


100 




1 


X 


10" 


27 


DQ Her 


PC 1015+14 


5940 


66 


160 


8 


X 


10" 


28 


IWD 


VV Pup 


6024 


145 


32 


3 


X 


10" 


28 


AM Her 


V834 Cen 


6090 


86 


22 


6 


X 


10" 


28 


AM Her 


MR Ser 


6816 


112 


25 


3 


X 


10" 


28 


AM Her 


BL Hyi 


6816 


128 


30 


3 


X 


10" 


28 


AM Her 


ST LMi 


6834 


128 


16 


3 


X 


10" 


28 


AM Her 


EXO 023432-5232.3 


6876 


500 




8 


X 


10" 


29 


AM Her 


TW Pic 


7188 


500 




7 


X 


10" 


29 


DQ Her 


EXO 033319-2554.2 


7590 


250 


56 


1 


X 


10" 


28 


AM Her 


Feige 7 


7920 


49 


35 


6 


X 


10" 


28 


IWD 


AM Her 


11136 


75 


13 


2 


X 


10" 


28 


AM Her 


V1500 Cyg 


12060 


1200 




1 


X 


10" 


29 


AM Her 


PG 1031+234 


12240 


142 


1000 


9 


X 


10" 


-29 


IWD 
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one-thousandth of the gas pressure, i.e. A4/\W\ < 3 x 10"''. If the magnetic white dwarfs discussed here had 
solar-like progenitors, and their fields are relic, the expected signal would be at least a factor of ten lower 
than depicted in Figure 

However, magnetic white dwarfs may have descended from a distinct population of main sequence stars, 



possibly magnetically active Ap stars {e.g. Angel, Borra & Landstrcet 1981) which have surface fields of 
10'^ — 10''' G and may harbor stronger internal fields than the sun. Dzicmbowski fc Goode (1985) argue 



that the frequency splittings of the modes rapidly oscillating Ap stars are proportional to the ratio of the 
magnetic to the gas pressure averaged over the excited region of the star. 

Both the magnetic field and the rotation of the star lift the degeneracy of modes with the same 
value of I but different values of m. By comparing the amplitude of each component of the perturbed 



modes and assuming a value of the rotational splitting constant ( ^jcdoux 1951| ), [Matthews (1991)| find 
(-PmagZ-Pgas) ~ 10"'^. Howcvcr, Dzicmbowski & Goode (1996) argue that the frequency splitting for the 



modes typically measured depends on the strength of the magnetic field in the outermost two percent 

radially of the star. To measure the global value of (Pmag/-fgas) and thereby obtain an estimate of 

one needs to study an ensemble of modes which together excite the star globally. Thus a similar difficulty 

arises in estimating the internal fields of the progenitors of magnetic white dwarfs as did for the white dwarfs 

themselves. 

Whatever has imparted an unusually strong magnetic field on the surfaces of the magnetic white dwarfs 
is likely to have boosted their internal fields as well. The measurement of gravitational radiation from rapidly 
rotating white dwarfs is a straightforward way to constrain the internal fields of these stars and provide a 
clue as to how their strong fields form. 

The author would like to thank Yu. Levin and E. S. Phinney for useful discussions and J. M. McDonald 
and L. Hernquist for providing comments on the manuscript. The author received support from a Lee A. 
DuBridge postdoctoral fellowship. 
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Fig. 1. — The reduced quadrupole moment Izz as a function of central density (pc) and magnetic field 
structure, p is the mean density of the white dwarf, W is the total gravitational potential energy of the 
star and Ai is the total energy of the internal magnetic field. T is the toroidal component of the magnetic 
energy. Each point represents an eigenstate of the coupled magnetic and stellar structure equations. 
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Fig. 2. — The estimated strain amplitude (h) of gravitational radiation from nearby white dwarfs. The solid 
squares are DQ Her variables, the empty squares are AM Her variables, the triangles are isolated magnetic 
white dwarfs and the asterisk is the dwarf nova/DQ Her WZ Sge. The second strongest source is AE Aqr 
which has a legend that depicts the scaling with the mass of the white dwarf on the right and with the 
root-mean-squared magnetic field strength on the left. The errorbars denote the theoretical uncertainty in 
h assuming that mass, radius, total magnetic energy and distance of the white dwarf are known precisely. 
The lower curve shows the sky-averaged noise level anticipated for a one-year LISA observation. The upper 
curve denotes a signal-to-noise ratio of five for a one-year observation. The dashed curve gives the binary 



confusion noise estimate (Bender ct al. 1998) 



